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Abstract. This paper presents the basic properties of chromatic derivatives 
and chromatic expansions and provides an appropriate motivation for intro- 
ducing these notions. Chromatic derivatives are special, numerically robust 
linear differential operators which correspond to certain families of orthogo- 
nal polynomials. Chromatic expansions are series of the corresponding special 
functions, which possess the best features of both the Taylor and the Shan- 
non expansions. This makes chromatic derivatives and chromatic expansions 
applicable in fields involving empirically sampled data, such as digital signal 
and image processing. 



1. Extended Abstract 

Let BL(7r) be the space of continuous L 2 functions with the Fourier transform 
supported within [— tt, tt] (i.e., the space of 7r band limited signals of finite energy), 
and let P£(cj) be obtained by normalizing and scaling the Legendre polynomials, 
so that 

^ J" P„» (w)dw = S(m - n). 

We consider linear differential operators JC n = (—i) n P£ (i gj); for such operators 
and every / e BL(n), 

fC n [m) = ^ f_ i B i^H/He^dc. 

We show that for / e BL(-7r) the values of K, n [f](t) can be obtained in a nu- 
merically accurate and noise robust way from samples of f(t), even for differential 
operators JC n of high order. 

Operators K7 1 have the following remarkable properties, relevant for applications 
in digital signal processing. 

Proposition 1.1. Let f : R — * R be a restriction of any entire function; then the 
following are equivalent: 

(a) £~o£ n [/](0) 2 <oc; 

(b) for allt^Jl the sum X^^Lg /C n [f](t) 2 converges, and its values are independent 
of t e R; 

(c) /gBL(tt). 
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Moreover, the following Proposition provides local representation of the usual 
norm, the scalar product and the convolution in BL(7r). 

Proposition 1.2. For all f, g e BL(7r) the following sums do not depend on t e R, 
and 



poo 

£/C"[/](t) 2 = / f(x) 2 dx; 

00 />oo 

Y,K n [f]{t)1C n [g\{t) = / f(x)g(x)dx; 

°c I-OO 

J2^ n lf]mt[9(u-t)} = / f(x)g(u-x)dx. 

The following proposition provides a form of Taylor's theorem, with the dif- 
ferential operators JC n replacing the derivatives and the spherical Besscl functions 
replacing the monomials. 

Proposition 1.3. Let j n be the spherical Bessel functions of the first kind; then: 

(1) for every entire function f and for all z£C, 

00 00 
f(z) = £(-l)"/C"[/](0)/C"[joM] - £/C n [/](0) v^T + T jnM; 

n=0 n=0 

(2) if f G BL(7r) 7 then the series converges uniformly on R and in L 2 . 

We give analogues of the above theorems for very general families of orthogonal 
polynomials. We also introduce some nonseparable inner product spaces. In one of 
them, related to the Hermite polynomials, functions /^(t) = sinwi for all w > have 
finite positive norms and for every two distinct values lo\ ^ loi the corresponding 
functions fu !l (t) = sinwii and fui 2 (t) = sinc^ are mutually orthogonal. Related to 
the properties of such spaces, we also make the following conjecture for families of 
orthonormal polynomials. 

Conjecture 1.4. Let P n (ui) be a family of symmetric positive definite orthonormal 
polynomials corresponding to a moment distribution function a(ui), 

/ -P„M -PmM da(uj) = S(m - n), 

J — oc 

and let 7„ > be the recursion coefficients in the corresponding three term recur- 
rence relation for such orth onormal polynomials, i.e., such that 

P n+1 (uj) = — P n (uj) - ^£zl P n -i(w). 
If In satisfy < lim — < 00 for some < p < 1, then 



n — 1 

< lim -j— V P k {iof < 00 

n—>oo n L ~P — ' 



n 1 

k=0 



for all u> in the support sp(a) of a(u>). 

Numerical tests with j n — n p for many p e [0, 1) indicate that the conjecture is 
true. 
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2. Motivation 

Signal processing mostly deals with the signals which can be represented by con- 
tinuous L 2 functions whose Fourier transform is supported within [— it, it]; these 
functions form the space BL(7r) of it band limited signals of finite energy. Foun- 
dations of classical digital signal processing rest on the Whittaker-Kotel'nikov- 
Nyquist-Shannon Sampling Theorem (for brevity the Shannon Theorem): every 
signal / € BL(-7r) can be represented using its samples at integers and the cardinal 
sine function sinci = short /irt, as 

oo 

(1) }{t)= J2 /(n) sine (t-n). 

n— — oo 

Such signal representation is of global nature, because it involves samples of the 
signal at integers of arbitrarily large absolute value. In fact, since for a fixed t the 
values of sine (t — n) decrease slowly as |n| grows, the truncations of the above series 
do not provide satisfactory local signal approximations. 

On the other hand, since every signal / G BL(7r) is a restriction to R of an entire 
function, it can also be represented by the Taylor series, 

(2) /w = £/ (n) (o)- 

n=0 

Such a series converges uniformly on every finite interval, and its truncations provide 
good local signal approximations. Since the values of the derivatives /^(O) are 
determined by the values of the signal in an arbitrarily small neighborhood of zero, 
the Taylor expansion is of local nature. In this sense, the Shannon and the Taylor 
expansions are complementary. 

However, unlike the Shannon expansion, the Taylor expansion has found very 
limited use in signal processing, due to several problems associated with its appli- 
cation to empirically sampled signals. 

(I) Numerical evaluation of higher order derivatives of a function given by its 
samples is very noise sensitive. In general, one is cautioned against numerical 
differentiation: 

". . . numerical differentiation should be avoided whenever possible, 
particularly when the data are empirical and subject to appreciable 
errors of observation" [10]. 
(II) The Taylor expansion of a signal / 6 BL(7r) converges non-uniformly on R; 

its truncations have rapid error accumulation when moving away from the 

center of expansion and are unbounded. 
(Ill) Since the Shannon expansion of a signal / € BL(7r) converges to / in BL(7r), 

the action of a continuous linear shift invariant operator (in signal processing 

terminology, a filter) A can be expressed using samples of / and the impulse 

response A [sine] of A: 

oo 

(3) A[f](t)= J2 f(n)A[smc](t-n). 

n— — oo 

In contrast, the polynomials obtained by truncating the Taylor series do not 
belong to BL(-7r) and nothing similar to ((3]) is true of the Taylor expansion. 
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Chromatic derivatives were introduced in [11] to overcome problem (I) above; the 
chromatic approximations were introduced in [14] to obtain local approximations 
of band-limited signals which do not suffer from problems (II) and (III) . 

2.1. Numerical differentiation of band limited signals. To understand the 
problem of numerical differentiation of band-limited signals, we consider an arbi- 
trary / 6 BL(7r) and its Fourier transform f(u>); then 

f (n \t) = i- J\u>rf{u)e^dw. 

Figure[T] (left) shows, for n — 15 to n — 18, the plots of (u)/%) n , which are, save a 
factor of i™, the symbols, or, in signal processing terminology, the transfer functions 
of the normalized derivatives l/V" d n /dt n . These plots reveal why there can be no 
practical method for any reasonable approximation of derivatives of higher orders. 
Multiplication of the Fourier transform of a signal by the transfer function of a 
normalized derivative of higher order obliterates the Fourier transform of the signal, 
leaving only its edges, which in practice contain mostly noise. Moreover, the graphs 
of the transfer functions of the normalized derivatives of high orders and of the 
same parity cluster so tightly together that they are essentially indistinguishable; 
see Figure ED (left) Q 

However, contrary to a common belief, these facts do not preclude numerical 
evaluation of all differential operators of higher orders, but only indicate that, from 
a numerical perspective, the set of the derivatives {/,/',/",...} is a very poor base 
of the vector space of linear differential operators with real coefficients. We now 
show how to obtain a base for this space consisting of numerically robust linear 
differential operators. 

2.2. Chromatic derivatives. Let polynomials i^(w) be obtained by normalizing 
and scaling the Legendre polynomials, so that 

i- £ P,» P* (w)(L; = S(m - n). 
We define operator polynomial^ 

JC? = H)^ (i • 

Since polynomials P£(u>) contain only powers of the same parity as n, operators 
K7 1 have real coefficients, and it is easy to verify that 

/C"[e iwt ] = i n P^{w)e lUt . 

Consequently, for / £ BL(7r), 

In particular, one can show that 

(4) /C n [sinc](f) = (-!)" V2^+lj„(7rt), 



If the derivatives are not normalized, their values can be very large and are again determined 
essentially by the noise present at the edge of the bandwidth of the signal. 

2 Thus, obtaining /C™ involves replacing ia k in P£(uj) with \ k d fe /dt fc for all k < n. If is applied 
to a function of a single variable, we drop index t in /C™. 




where j n (x) is the spherical Bessel function of the first kind of order n. Figure [T] 
(right) shows the plots of P£(u>), for n — 15 to n = 18, which are the transfer 
functions (again save a factor of i") of the corresponding operators JC n . Unlike the 
transfer functions of the (normalized) derivatives I/ti 71 d"/di ra , the transfer func- 
tions of the chromatic derivatives K. n form a family of well separated, interleaved 
and increasingly refined comb filters. Instead of obliterating, such operators encode 
the features of the Fourier transform of the signal (in signal processing terminol- 
ogy, the spectral features of the signal) . For this reason, we call operators K n the 
chromatic derivatives associated with the Legendre polynomials. 

Chromatic derivatives can be accurately and robustly evaluated from samples 
of the signal taken at a small multiple of the usual Nyquist rate. Figure [2] (left) 
shows the plots of the transfer function of a transversal filter given by 7^5 [/] (t) = 
Sfcl-64 Cfe (g ra y)j used to approximate the chromatic derivative /C 15 [/](i), 

and the transfer function of /C 15 (black). The coefficients Ck of the filter were ob- 
tained using the Remez exchange method [16], and satisfy \ck\ < 0.2, (—64 < k < 
64). The filter has 129 taps, spaced two taps per Nyquist rate interval, i.e., at a 
distance of 1/2. Thus, the transfer function of the corresponding ideal filter /C 15 
is P 1 i 5 (2o;) for |u;| < tt/2, and zero outside this interval. The pass-band of the ac- 
tual transversal filter is 90% of the bandwidth [—tt/2, tt/2]. Outside the transition 
regions [— ll7r/20, — 97r/20] and [97r/20, ll7r/20] the error of approximation is less 
than 1.3 x 10~ 4 . 

Implementations of filters for operators JC n of orders < n < 24 have been tested 
in practice and proved to be both accurate and noise robust, as expected from the 
above considerations. 

For comparison, Figure[2] (right) shows the transfer function of a transversal filter 
obtained by the same procedure and with the same bandwidth constraints, which 
approximates the (normalized) "standard" derivative (2/-7r) 15 d 15 /di 15 (gray) and 
the transfer function of the ideal filter (black). The figure clearly indicates that 
such a transversal filter is of no practical use. 

Note that (JD) and (@} imply that 

(5) 

OO OO 

£ k im)= E f(n)IC k [smc}(t-n)= ]T f(n) (-l) k V2kTT j fc (ir(t - n)). 
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FIGURE 2. Transfer functions of JC 15 (left, black) and d 15 /dt 15 
(right, black) and of their transversal filter approximations (gray). 
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Figure 3. left: Oscillatory behavior of sine (t) (black), and 
/C 15 [sinc](i) (gray); right: A signal / E BL(7r) (gray) and its 
chromatic and Taylor approximations (black, dashed) 



However, in practice, the values of K, k [f](t), especially for larger values of k, cannot 
be obtained from the Nyquist rate samples using truncations of ([5]). This is due to 
the fact that functions /C fc [sinc](t — n) decay very slowly as \n\ grows; see Figure [3] 
(left). Thus, to achieve any accuracy, such a truncation would need to contain 
an extremely large number of terms. On the other hand, this also means that 
signal information present in the values of the chromatic derivatives of a signal 
obtained by sampling an appropriate filterbank at an instant t is not redundant 
with information present in the Nyquist rate samples of the signal in any reasonably 
sized window around t, which is a fact suggesting that chromatic derivatives could 
enhance standard signal processing methods operating on Nyquist rate samples. 



2.3. Chromatic expansions. The above shows that numerical evaluation of the 
chromatic derivatives associated with the Legendre polynomials does not suffer 
problems which precludes numerical evaluation of the "standard" derivatives of 
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higher orders. On the other hand, the chromatic expansions, defined in Proposi- 
tion [2J] below, were conceived as a solution to problems associated with the use of 
the Taylor expansion 

Proposition 2.1. Let JC n be the chromatic derivatives associated with the Legendre 
polynomials, let j n be the spherical Bessel function of the first kind of order n, and 
let f be an arbitrary entire function; then for all z,u€ C, 

oo 

(6) f(z) = K n [f]{u)K^mc{z-u)\ 

oo 

(7) = ^(-l)"^"[/](«)^"[sinc](z- M ) 

71=0 



(8) = ^2K n [f](u)V2^+lj n (^-u)) 

n=a 



If f G BL(tt), then the series converges uniformly on R and in the space BL(7r). 

The series in ([6]) is called the chromatic expansion of f associated with the Le- 
gendre polynomials; a truncation of this series is called a chromatic approximation 
of /. As the Taylor approximation, a chromatic approximation is also a local ap- 
proximation; its coefficients are the values of differential operators IC m [f}(u) at a 
single instant u, and for all k < n, 



jfc 

/ (fc) («) = TT 



JC m [f](u) ACS 1 [sine (t-«)] 



Figure[3] (right) compares the behavior of the chromatic approximation (black) of 
a signal / 6 BL(7r) (gray) with the behavior of its Taylor approximation (dashed). 
Both approximations are of order 16. The signal f(t) is defined using the Shannon 
expansion, with samples {f(n) : \f(n)\ < 1, —32 < n < 32} which were randomly 
generated. The plot reveals that, when approximating a signal / € BL(7r), a 
chromatic approximation has a much gentler error accumulation when moving away 
from the point of expansion than the Taylor approximation of the same order. 

Unlike the monomials which appear in the Taylor formula, functions K, n [sine] (t) = 
(-1)" ^2n + 1 j„(7rf) belong to BL(rr) and satisfy |£"[sinc](t)| < 1 for all t S R. 
Consequently, the chromatic approximations also belong to BL(7r) and are bounded 
on R. 

Since by Proposition 12.11 the chromatic approximation of a signal / G BL(7r) 
converges to / in BL(tt), if A is a filter, then A commutes with the differential 
operators JC n and thus for every / £ BL(tt), 

oo 

(9) A[f](t) = ^(-l) n JC n [f](u) K n [A[smc]}(t - u). 

71 = 

A comparison of ([9]) with ([3]) provides further evidence that, while local just like 
the Taylor expansion, the chromatic expansion associated with the Legendre poly- 
nomials possesses the features that make the Shannon expansion so useful in sig- 
nal processing. This, together with numerical robustness of chromatic derivatives, 



^Propositions stated in this Introduction are special cases of general propositions proved in sub- 
sequent sections. 
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makes chromatic approximations applicable in fields involving empirically sampled 
data, such as digital signal and image processing. 

2.4. A local definition of the scalar product in BL(7r). Proposition ^. 21 below 
demonstrates another remarkable property of the chromatic derivatives associated 
with the Legendre polynomials. 

Proposition 2.2. Let f : R — > R be a restriction of an arbitrary entire function; 
then the following are equivalent: 

(b) for allt(ER the sum X)n=o ^-"[/] W 2 converges, and its values are independent 
of t G R; 

(c) /gBL(tt). 

The next proposition is relevant for signal processing because it provides local 
representations of the usual norm, the scalar product and the convolution in BL(7r), 
respectively, which are defined globally, as improper integrals. 

Proposition 2.3. Let)C n be the chromatic derivatives associated with the (rescaled 
and normalized) Legendre polynomials, and f,g£ BL(tt). Then the following sums 
do not depend on t G R and satisfy 

(10) E^™*) 2 = / Kxfdx; 

n=0 J ~°° 
°° y* oo 

(11) £* n [/](t)* T Mt) = / f(x)g(x)dx; 

°° rOO 

(12) Y, Kn ^ K t^-t)]= f{x)g{u-x)dx. 

n=0 J -°° 

2.5. We finish this introduction by pointing to a close relationship between the 
Shannon expansion and the chromatic expansion associated with the Legendre poly- 
nomials. Firstly, by (J7J), 



(13) f(n) = £ ^ fe [/](0) (-l) fc /C fc [sinc](n). 

fc=0 

Since ZC n [sinc](£) is an even function for even n and odd for odd n, ([5|) implies 

oo 

(14) /C fc [/](0)= J2 f(n)(-l) k K k [ S inc](n). 



TI— — OC 



Equations (fT3")l and (jT^J) show that the coefficients of the Shannon expansion of a 
signal - the samples f(n), and the coefficients of the chromatic expansion of the 
signal - the simultaneous samples of the chromatic derivatives /C n [/](0), are related 
by an orthonormal operator defined by the infinite matrix 



[(-l) fe /C fe [sine] (n) : keW,ne'l]= V2fc + 1 j fc (7m) : fcei,!i£ 

Secondly, let S u [f(u)] = f(u + 1) be the unit shift operator in the variable u 
(/ might have other parameters). The Shannon expansion for the set of sampling 
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points {u + n : n G Z} can be written in a form analogous to the chromatic 
expansion, using operator polynomials <S™ = S u ° ■ ■ ■ ° S u , as 

oo 

(15) f(t) = ^/( U + n)sinc(i-( U + n)) 

71=0 
OO 

(16) = £ £?[/](«) S>mc(i-«)]; 

n=0 

compare now (|16|) with ((6]). Note that the family of operator polynomials {5™}„ e ^ 
is also an orthonormal system, in the sense that their corresponding transfer func- 
tions {e inw } n6Z form an orthonormal system in L 2 [— 7r,7r]. Moreover, the transfer 
functions of the families of operators {/C"}„ e]N and {6>™} nez , where K. n are the 
chromatic derivatives associated with the Legendre polynomials, are orthogonal on 
[— 7r, 7r] with respect to the same, constant weight w(u>) — l/(27r). 

In this paper we consider chromatic derivatives and chromatic expansions which 
correspond to some very general families of orthogonal polynomials, and_prove gen- 
eralizations of the above propositions, extending our previous work [12] □ However, 
having in mind the form of expansions ((6]) and (I16|) . one can ask a more general 
(and somewhat vague) question. 

Question 1. What are the operators A for which there exists a family of operator 
polynomials {P n (A)}, orthogonal under a suitably defined notion of orthogonality, 
such that for an associated function m A (t), 

f(t) = ^P n (A)[/](u) P n (A)[m A (t-u)} 

n 

for all functions from a corresponding (and significant) class C A ? 

3. Basic Notions 

3.1. Families of orthogonal polynomials. Let M. : Vu — ► 1R be a linear func- 
tional on the vector space of real polynomials in the variable w and let fi n — 
M(u> n ). Such M. is a moment functional, fi n is the moment of M. of order n and 
the Hankel determinant of order n is given by 





Mo ■ 




A n = 


Mi ■ 


■ Mn+1 




[in ■ 


M2n 



The moment functionals M. which we consider are assumed to be: 

(i) positive definite, i.e., A„ > for all n; such functionals also satisfy fi2n > 0; 

(ii) symmetric, i.e., /^n+i = for all n; 

(iii) normalized, so that = /j,q = 1. 

For functionals M. which satisfy the above three conditions there exists a family 
{P^(uj)}n£N of polynomials with real coefficients, such that 

Chromatic expansions corresponding to general families of orthogonal polynomials were first 
considered in [5]. However, Proposition 1 there is false; its attempted proof relies on an incorrect 
use of the Paley- Wiener Theorem. In fact, function F g defined there need not be extendable to 
an entire function, as it can be shown using Example 4 in Section[5]of the present paper. 
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(a) {P T ^ 4 (o;)} rie ] N is an orthonormal system with respect to M, i.e., for all m, n, 

M{P£{u>)P?{u>)) = 8{m-n); 

(b) each polynomial P™{oj) contains only powers of u of the same parity as n; 

(c) F-H = l. 

A family of polynomials is the family of orthonormal polynomials corresponding 
to a symmetric positive definite moment functional M. just in case there exists a 
sequence of reals 7„ > such that for all n > 0, 



(17) C + lH = -Pn{u) ~ ^P,f-lM- 

7™ ]n 

If we set 7_i = 1 and P^j(w) = 0, then (fTT)) holds for n = as well. 

We will make use of the Christoffel-Darboux equality for orthogonal polynomials, 



(18) ( W - a) £ P£*{")Pk (*) = 7n(i^iHiT to ~ i^-iOWM), 

fe=0 

and of its consequences obtained by setting a = —lo in (fT8|) to get 

(19) u, J^+iM 2 - E ) = 72„+i Pg+M ^M, 

Vfe=o fc=o / 

and by by letting <r — > a; in (|18p to get 

n 

(20) £ p fc " m 2 = 7 „ TOxM'p^M - /r+iM^f M')- 

fc=0 

For every positive definite moment functional M. there exists a non-decreasing 
bounded function a(u>), called a moment distribution function, such that for the 
associated Stieltjes integral we have 



(21) / uj n da{uj) = fi n 

J — OO 

and such that for the corresponding family of polynomials {P r f 1 (a;)} ne ]>j 

(22) f" P^{u) da(w) - 5(m - »). 



We denote by L 2 a ^ the Hilbert space of functions <p : R — > C for which the 
Lebesgue-Stieltjes integral |<,£>(ci;)| 2 da(w) is finite, with the scalar product de- 
fined by (a, (3) a (ui) — ct(uj) (3(u>) da(u>), and with the corresponding norm de- 
noted by |M| o(w) . 

We define a function m : R — > R as 

/OO 
e iuJt da(uj). 
-OO 

Since 

/•oo / /» oo ^oo \ 1/2 

/ \(iu) n e iwt \da(uj) < ( / uj 2n da(u) do(w) ) = x/T^T < oo, 
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we can differentiate (|23|) under the integral sign any number of times, and obtain 
that for all sn, 



(24) 
(25) 



m ( 2 ")(0) = (-1)"M2„; 

m (2n+l)( ) = 0. 



3.2. The chromatic derivatives. Given a moment functional M. satisfying con- 
ditions (JTJ) - jm{ above, we associate with Ai a family of linear differential operators 
{/C"}^ 6]N defined by the operator polynomial 



K2 



and call them the chromatic derivatives associated with Ai. Since Ai is symmetric, 
such operators have real coefficients and satisfy the recurrence 



(26) 



K. n+1 = —{DoK, n ) + — IC n -\ 

In In 



with the same coefficients 7„ > as in (fT71) . Thus, 
(27) /C?[e iwt ] = i n P^{w)e iu,t 

and 



(28) 



K, n [m]{t) 



i"P^(w)e i ^da(a;). 



The basic properties of orthogonal polynomials imply that for all to, n, 

(29) (tC n o /C m )[m](0) = {-l) n S(m - n), 
and, if to < n or if to — n is odd, then 

(30) (£» m o/C")[m](0) = 0. 

The following Lemma corresponds to the Christoffel-Darboux equality for or- 
thogonal polynomials and has a similar proof which uses (|26|) to represent the left 
hand side of (l3T|) as a telescoping sum. 

Lemma 3.1 ([12]). Let {JC n }n% K be the family of chromatic derivatives associated 
with a moment functional Ai, and let f,g G C°° ; then 



(31) 



D 



.m=0 



7 n(/C n+i [/] AC" [£]+£"[/]£ 



n+lr 



5 Thus, to obtain K. n , one replaces u) k in P**(o>) by i fc D k , where Dj [/] = 4rnf(t). We use the 
square brackets to indicate the arguments of operators acting on various function spaces. If A 
is a linear differential operator, and if a function f(t,w) has parameters w, we write At[/] to 
distinguish the variable t of differentiation; if f(t) contains only variable t, we write A[f(t)] for 
A t [/(*)] and D k [/(*)] for D k [/(*)]■ 
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3.3. Chromatic expansions. Let / be infinitely differentiable at a real or complex 
u; the formal series 

oo 

(32) CET\f,u](t) = E^*[/](«)^[n»(t-«)] 

oo 

= ^(-l) fc /C fc [/]( U )/C fe [m](i- U ) 

k=0 

is called the chromatic expansion of / associated with M, centered at u, and 

n 

CA M [/, n, u] (t) = £ (- 1 ) fc /C fc [/] ( M )/C fc [m] (i - u) 

fc=0 

is the chromatic approximation of / of order n. 

From (|29|) it follows that the chromatic approximation CA^ [/, n, u](t) of order 
n of /(t) for all in < n satisfies 

n 

K7*[CA"[/,n, «](*)] | t=u = ^(-l)V[/]H(ro^)[m](0) = /C m [/]( U ). 

fc=0 

Since /C™ 1 is a linear combination of derivatives D k for k < m, also f^ m ^(u) = 
Dt*[CA M [/, n, for all m < n. In this sense, just like the Taylor approxi- 

mation, a chromatic approximation is a local approximation. Thus, for all m < n, 

n 

(33) (u) = £>r [CA" [/, n, «](*)] | t=i = ^(-l) fc /C fc [/]( M ) (D m o JC k ) [m] (0) . 

fc=0 

Similarly, since [E*=o f W " «) fc A ! ] | t=u = / (m) («) for 

to < n, we also 

have 



(34) K, m [f]{u) = K,\ 



J2f (k) (u)(t-u) k /k\ 



.k=0 J t=u fe=0 



= ]T/«( U )/C m [t*/fc!] (0). 



Equations (|33|) and (|34[) for m = rt relate the standard and the chromatic bases 
of the vector space space of linear differential operators, 



(35) D n = ^(-l) fc OD"o/C fe )[m](0) JC k ; 

k=0 
n 

(36) K n = ^K n [t k /k\\{Q) D k . 

k=a 

Note that, since for j > k all powers of t in JC k ]t? are positive, we have 

(37) j>k K k [tii/jl] (0) = 0. 

4. Chromatic Moment Functionals 

4.1. We now introduce the broadest class of moment functionals which we study. 

Definition 4.1. Chromatic moment functionals are symmetric positive definite 
moment functionals for which the sequence {fin /n}neia * s bounded. 
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If M is chromatic, we set 



, QQ s y /Mn\ 1/n ,. fin" . 
(38) p = lira sup I — -J = e lim sup < oo. 

Lemma 4.2. Let M be a chromatic moment functional and p such that (|38[) holds. 
Then for every a such that < a < l/p the corresponding moment distribution 
a(u>) satisfies 



l/r, 



(39) / e a ^da{tu) < oo. 

J — OO 

Proof. For all b > 0, /* fc e Q Mda(tj) = £~ =0 /! b M n da(w). For even n we 

have J , oj n da(uj) < p n . For odd n we have < 1 + uj n+1 for all and thus 
J* , \uj\ n da(Lu) < J^ b da(uj) + J. uj n+1 da(uj) < 1 + ^n+i- Let /?„ = fi n if n is even, 
and f3 n — 1 + p n +i if is odd. Then also e a l aj 'da(a;) < arl /W n " Since 

limsup n ^ 00 (/3„/n!) 1 /" = p and < a < l/p, the last sum converges to a finite 
limit. □ 

On the other hand, the proof of Theorem 5.2 in §11.5 of [7] shows that if (|39|) 
holds for some a > 0, then (|38|) also holds for some p < 1/a. Thus, we get the 
following Corollary. 

Corollary 4.3. A symmetric positive definite moment functional is chromatic just 
in case for some a > the corresponding moment distribution function a(u>) satis- 
fies ([39]) . 



Note: For the remaining part of this section we assume that M. is a chromatic 
moment functional. 

For every a > 0, we let S(a) = {zeC : |Im(z)| < a}. The following Corollary 
directly follows from Lemma [4~2l 

Corollary 4.4. IfueS(^), then e iuw G L 2 , y 

We now extend function m(t) given by (f2"3"| from R to the complex strip S(^). 
Proposition 4.5. Let for z 6 S(-), 

/oo 
e iuJz da(uj). 
-oo 

Then m(z) is analytic on the strip S(j^). 

Proof. Fix n and let z — x + iy with \y\ < l/p; then for every b > 0, 

|(iw) n e i<J *|da(a;) < /Vr el ^'da(^) = V^T / M" +fc daM. 

As in the proof of Lemma 14.21 we let (3k = p n +k for even values of n + k, and 
[3k = 1 + A*n+fc+i for odd values of n + k; then the above inequality implies 

00 I I k R 

\(iu) n ^"\6a{u) <^ 



^ fc! 

fe=0 



and it is easy to see that for every fixed n, limsup fc ^ 00 (/3fc/fc!) 1 /''" = P- ^ 
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Proposition 4.6. Let (p(u>) <E L 2 ^; we can define a corresponding function f v 



(41) U(z)= / ip(uj)e MZ da(w). 

J — OO 

Such f<p{z) is analytic on S(j^) and for all n and z £ S(^-) } 

(42) 1C n [U]{z)= \ n P?{w)y{u)4»*&a{u). 



Proof. Let z — x + iy, with \y\ < 1/(2/9). For every n and b > we have 

\(iu) n (p(u)fP"\da(u) 
,6 

< / \u\ n \(p{u)\e^da(u) 

J-b 



fc=0 ' 

Qg I ife / fb r-b \ 



Thus, also 



/ |(ic)>Me^|daH < ||^M|| aM £ TTV^ 

J —OG i n 



2k- 

/,'! ' 

fc=0 



The claim now follows from the fact that limsupj,^,^ 2 ^//i2n+2fc/ = 2p for every 
fixed n. □ 

Lemma 4.7. //.M is chromatic, then {P^i^^neiN * s ° complete system in L^, u y 

Proof. Follows from a theorem of Riesz (see, for example, Theorem 5.1 in §11.5 
of [7]) which asserts that if liminf„_ ) . 00 {^i n /n\) 1 ^ n < oo, then {P^(u;)} n€ jN is a 
complete system in L^^H 

□ 

Proposition 4.8. Let ip(oj) <E -^o(w)' */ / or some ^ced u € S(jfi) the function 
(p(u>)e iuu also belongs to L 2 a ^, then in L 2 a ^ we have 



(43) = ^(-ij-rt/jt.) 

n=0 

and for f v {z) given by (|41j) we /iai/e 

oo 

(44) i^iuku)] 2 = hw^Wlw < oo. 



''Note that we need the stronger condition limsup n _ (00 (/im/n!) 1 '" < oo to insure that function 
m(z) denned by l|40jl is analytic on a strip (Proposition 14. 5) . 
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Proof. By Proposition 14.61 if u 6 S(j-), then equation ([4"2f holds for the corre- 



sponding f v given by ([4"Tjl . If also (,c(w)e iwu € F 2 ( w ), then (|4^|) implies that 



(45) (^He-",P r f H) aM = (-irXT^Ku). 

Since {P^ l (£j)}„g ][ f is a complete orthonormal system in ^/ ■., (j45|) implies (|43|) . 
and Parseval's Theorem implies (|44|) . □ 



Corollary 4.9. For every (p(u>) € ^o(a>) arl ^ every u£R, equality f|43[) ZioWs and 



(46) 



£|/C"[/ V ](u)| 2 = UpHIG 



n=0 



Taws, i/ie sum X^Lo l^™[/v]( u )| 2 * s independent of u £ R. 
Proo/. If it S R, then ip{u})e^ u e L 2 (w) and II^H^"^^) 
Corollary 4.10. Lei e > 0; then for all u e S(j- - e) 



IvHIl* 



□ 



(47) 



^|/C"[m]( U )| 2 < 



n=0 



e (^)M 



a(aj) 



< OO. 



Proof. Corollary 14.41 implies that we can apply Proposition 14.81 with ip(u>) = 1, in 
which case f v (z) = m(z), and, using Lemma 14.21 obtain 



J2\£ n lrn}(u)\ 2 = \\e^f a 



n=0 



< 



JIm(?i) oj 



a(uj) 



< 



Definition 4.11. L 2 , is i/ie vector space of functions f : S(j^) 
analytic on S(^) and satisfy X^^Lo I ^"[/l (0)1 2 < 00 • 
Proposition 4.12. The mapping 

OO 

wM = £(-i) n r i [/](o)i^H 



< OO. 

o(w) 

□ 

— ► C which are 



(48) 



/(*) 



71=0 



is an isomorphism between the vector spaces LjJ, and i 2 ^, and its inverse is given 



by 65). 

Proo/. Let / G LjJ,; since 1^-™ I/K^) I 2 < 00 j tne function </5/(w) defined by 



48]) belongs to F 2 ^)- By Proposition 14.61 f Vf defined from <pj by (jITj) is analytic 
on and by Proposition [U it satisfies <p f (w) = E^=o(-i) n ^[/*>/]( ) P rf M- 

By the uniqueness of the Fourier expansion of <Pf(ui) with respect to the system 
{P£»}„ e]N we have /C™[/](0) = /C n [/ v/ ](0) for all n. Thus, /(«) = f Vf (z) for all 

*eS(£). □ 

Proposition 14. 121 and Corollary 14.91 imply the following Corollary. 

Corollary 4.13. For all f 6 \? M and alltER the sum X^lo l^-"[/](^)| 2 converges 
and is independent oft. 

Definition 4.14. For every f(z) G L~, we call the corresponding tpf(oj) G 
given by equation (|48[) £/ie M- Fourier- Stieltjes transform of f(z) and denote it by 
^ [/]("). 
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Assume that a(u>) is absolutely continuous; then a'(w) = w(w) almost everywhere 
for some non- negative weight function w(cj). Then (|4T|) implies 

/•OO 

f(z)= / r M [f](aj)j»*w(<j)du>. 



This implies the following Proposition. 

Proposition 4.15 ([12]). Assume that a function f{z) is analytic on the strip 
S'(y-) and that it has a Fourier transform f(uj) = f{t) e~ iwt dt such that 

f( z ) = h J^/M eiuJZ duJ f° r aU z e then f( z ) e L li if and onl v 

J^lo l/^)! 2 w (^) _1 < oo, in which case f(u>) = 27r J 7 ^ [/](cj) w(w). 

4.2. Uniform convergence of chromatic expansions. The Shannon expansion 
of an / 6 BL(7r) is obtained by representing its Fourier transform f(ui) as series 
of the trigonometric polynomials; similarly, the chromatic expansion of an / G 
Lj^ is obtained by representing J rM [f](u>) as a series of orthogonal polynomials 

{p r f MW 

Proposition 4.16. Assume f S Lj^; i/iew /or all m £ R anc! e > 0, £/ie chromatic 
series CE M [/, u](z) of f{z) converges to f(z) uniformly on the strip S(j^ — e). 

Proof. Assume it € R; by applying (|42[) to m(z — u) we get that for all z G S(-%-), 

/oo 71 
-°° fe=0 

Since / G L£, Proposition 14.121 implies T M [f](uj) G Corollary Ol and 

equation (|43|) imply that in i^/^s 

oo 

e 1 - = E(-i) fc £ fc [/](«)n M M- 

fc=0 

Thus, 

/OO 00 
E(-i)^ fc [/](")^(^)e iw( ^ ll) da(c). 
-oo 



Consequently, from (|49|) and {5 

|/(*)-CA"[/,n,u](«)| 



< 



E (-i) fe /C fe [/]( M )P fc -Me U 



fe=n+l 



< 



E (-i) n £"[/](«WM 
For z G — e) we have 



(2- t() 



da (a;) 



da(u>) 



1/2 



= iw(z— u) |2 



da(w) 



(51) \f(z) -CA M [f,n,u](z)\ < ( E l^[/](")| 2 / 

\fc=n+l ^ 



1/2 



e (i-2e)|^| 



da(w) 
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Consequently, Lemma l4~2l and Corollary 14. 131 imply that CE A4 [/,u](z) converges to 
f(z) uniformly on S(j- — e). □ 

Proposition 4.17. Space consists precisely of functions of the form f(z) = 
yy^Ln Qn/C" [wt| (z) where a = ((a„))„ e]N is a complex sequence in I 2 . 

Proof. Assume a £ I 2 ; by Proposition ^. 10[ for every e > 0, if z 6 S{^- — e) then 

oo / oo oo \ 

X> n /C»M(s)| < ^|a„| 2 ^|/C"[m](z)| 2 

n=fe \ n — /c n— /c / 

/ oo \ V2 



< 



which implies that the series converges absolutely and uniformly on S(j^ — e). 
Consequently, f(z) — o- n K- n [m\{z) is analytic on S(j-), and 

oo 

£" l [/](0) =^«„r^")[m](0) = (-l) m a m . 

n=Q 

Thus, En=o\£ n [/](«) I 2 = E,"oKI 2 < 00 and 80 /(*) e L*. Proposition OH 
provides the opposite direction. □ 

Note that Proposition 14. 1 71 implies that for every e > functions f(z) E L 2 ^ are 
bounded on the strip S(j- — e) because 



i/wi< £™(o)i 2 



1/2 






e (^-«0M 



4.3. A function space with a locally defined scalar product. 

Definition 4.18. is the space of functions f(t) : R i— * C obtained from functions 
in L 2 ^ by restricting their domain to R. 

Assume that f,gG Lw ; then (|4"5)) implies that for all m £ R, 



n=0 



Note that for all t, u e R, 

n 



ri=0 



1 fc=0 



By (031, the sum ££ =0 (- i)*£ fc [/] W converges in L 2 (tij) to T M 

Since ^ [#](w) e^C-") G i 2 (w) , and since 



^[ 5 ]( W )^[/]( W )e 1 - t da( W ) 



<ll^[/]HII BH ||^[ff]HH o(w) 
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we have that for all t, u £ R, 

oo /*°° 

£/C fe [/] (u) - U )] = / ^[ S ](w)^[/](w)e^ t da(w)<oo. 

I n — OO 



Proposition 4.19 ([12]). We con introduce locally defined scalar product, an as- 
sociated norm and a convolution of functions in Lj^ by the following sums which 
are independent of u £ R : 

OO 

(52) \\f\t = El^ n [/]^)I 2 = II^[/]HII" M ; 

oo 

(53) (f,g) M = "£K n [f](u)Kn[g](u) 

n=0 

= (^[/]M,^[ 5 ]H) aM ; 



(54) (f* M 9)(t) = J2 Kn lf^ K ul9(t-u)} 

T M [f]{uj) ^[ 5 ](a;)e^*da(a;). 



OO 



Letting g(t) = m(t) in $5U) , we get (/ * M m){t) = C&^[f,u]{t) = f{t) for all 
f(t) £ L^, while by setting u = 0, u = t and u = t/2 in (|54|) . we get the following 
lemma. 

Lemma 4.20 ( [12]). For every f,g£ and for every teR, 

OO OO 

J2(-l) k IC k [f}(t)IC k [g}(Q) = ^(-l) fe /C fc [/](0)/C fc [ 3 ]W 

k=0 k=0 

oo 

= ^(-l) fc /C fc [/](t/2)/C fc [.g](t/2). 

fc=0 

Since m(z) is analytic on S(-), so are /C"[m](z) for all n; thus, since by (j2"5|) 

Em,o( roK '1 ffl l(°) 2 = we have ^"MW G L £ for a11 n - Let u be a fixed real 
parameter; consider functions B™(t) = /C"[m(i — u)] = (— 1)™/C™[m](t — u) £ 1/^. 

Since 
(55) 

OO 

= £(/C t fc o ^)[m(i - «)](/C t fe o K%)[m(t - «)] 

oo 

= ^(- 1 ) m+n ( /cfe ° - u){K. k o /C m )[m](i - u) 

= 5(m — n), 

the family {_B™(t)}„ e]N is orthonormal in LjJ, and for all n £ IN and all t £ R, 

OO 

(56) ^2{JC k o JC n )[m}{tf = 1. 

fc=0 
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-2 
-'m i 



By ([H, for / € L 

(f,KZ[m(t-u)]) M = £/C fe [/](i)(/C t fc o/C£)[m(i- U )] 



fc=0 



£(-l) n /C fc [/]( U )(/C fc o/C")[m](0) 



k=0 

(57) = K n [f](u). 

Proposition 4.21 ([12]). The chromatic expansion CP3 M [/, u](t) of f(t) £ Lj^ is the 
Fourier series of f(t) with respect to the orthonormal system {/C"[m(i — u)]}neiN- 
The chromatic expansion converges to f(t) in L~^; thus, {/C™[m(i — u)]}„ 6 ]n is a 
complete orthonormal base of LjJ,. 

Proof. Since JC k [f{t) - CET 1 [/, n, u](t)] | t=u equals for jfe < n and equals /C fe [/](u) 
forfc>n, ||/-CEr[/,n )W ]|| JM =Er=n+i^[/]M 2 ^0- D 

Note that using (|56|) with 7i = we get 

oo 

|/(f)-CE M [/,n J «](t)| < £ |/C fc [/]( M )/C fc [m](t- U )| 

fc=n+l 

/ oo oo \ V 2 

< 53 K k [f]{uf £ K k lrn](t-uf\ 

\k=n+l fe=n+l / 

/ oo \V2 / V 1/2 

(58) = £ £*[/](«) a l-£/C fc [m](t-n) 2 . 

\fe=n+l / V fe=0 / 

Let 

/ n \ V2 

Bn(i)= l-£/C fc [m](t) 2 ; 
V fe=0 / 
then, using Lemma 13. 11 we have 

/ n \~V2 
K(t) = 7 « /C n+1 [m](t) /C n [m](i) 1 - £/C fe [m](i) 2 

V fe=0 / 

Since (D k o /C n )[m](0) = for all < k < n - 1, we get that E [k \t) = for all 
k < 2n + 1. Thus, E n (0) = and E n (t) is very flat around f = 0, as the following 
graph of £15 (t) shows, for the particular case of the chromatic derivatives associated 
with the Legendre polynomials. This explains why chromatic expansions provide 
excellent local approximations of signals / € BL(7r). 

4.4. Chromatic expansions and linear operators. Let A be a linear operator 
on LjJ, which is continuous with respect to the norm U/H^. If A is shift invariant, 
i.e., if for every fixed h, A[f(t + h)] = + h) for all / € L 2 ^, then A commutes 

with differentiation on and 

00 

A[f](t) = £ {-l) n K n [f]{u) JC n [A[m]](t - u) = (/ * M A[m])(t), 

n=0 
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Figure 4. Error bound £15 (i) for the chromatic approximation 
of order 15 associated with the Legendre polynomials. 



Consequently, the action of such A on any function in LjJ, is uniquely determined 
by A[m], which plays the role of the impulse response A[sinc] of a continuous 
time invariant linear system in the standard signal processing paradigm based on 
Shannon's expansion. 

Note that if A is a continuous linear operator A on such that (AoD n )[m] (t) = 
(D n o A)[m](t) for all n, then Lemma T4.20I implies that for every / G L^, 

oo 

A[f](t) = £(-l) n JC»[/](0)r*[A[m]](t) 

OO 

= £(-l)»;C»[A[m]](0) £»[/](*). 

n=0 

Since operators /C n [/](t) are shift invariant, such A must be also shift invariant. 

4.5. A geometric interpretation. For every particular value of t G IR the map- 
ping of into I 2 given by / h- ► f t — ((£"[/](£) ))neN is unitary isomorphism which 
maps the base of L^, consisting of vectors B k (t) = (—l) k K, k [m(t)], into vectors 
B k = (((-l) k (JC n o jc k )[m(t)])) n£K . Since the first sum in ((SSJ) is independent of t, 
we have {B k ,B r t n ) = S(m-k), and ([57]) implies (/ f ,B t fe ) = K k [f](0). Thus, since 
££o£*[/](0) a < oo, we have £~ ^ M(0) B k G ; 2 and 

OO OO 

£</ t ,B*)B* = £x:*[/](o)b* 

fc=0 fc=0 

oo 

= ((^/C fc [/](0)(-l) fe (/C"o/C fc )[m(t)])) neK . 

Since for / e Lj^ the chromatic series of / converges uniformly on R, we have 
K n [f](t) = EZoK k [f]m-l) k CC n -JC k )[m(t)]. Thus, 

oo oo 

J2(fuB k ) B k =J2K k [fW) B k = {{K, n [f]{t))) nen = f t . 

Thus, while the coordinates of ft = ((/C™ [f](t) )) n eiN m the usual base of I 2 vary 
with t, the coordinates of f t in the bases {-B^j-fce^ remain the same as t varies. 
We now show that {£?™}„ e]N is the moving frame of a helix H : R i— > i 2 - 
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Lemma 4.22 ([12]). Let / € I? M and let t € R vary; then /(*) = ((/C n [f](t))) neK 
is a continuous curve in I2. 



Proof. Let / e Ll 



then, since YlnLo^ 



[f]{t) 2 converges to a continuous (con- 



stant) function, by Dini's theorem, it converges uniformly on every finite interval 
/. Thus, the last two sums on the right side of inequality \\f(t) — f(t + h)\\ M < 

Eto(£ n [/](*) - + h )f +2En= N+ i ^ n \m? +2E: w K n [/](i + 

can be made arbitrarily small on I if N is sufficiently large. Since functions 
)C n [f](t) have continuous derivatives, they are uniformly continuous on /. Thus, 
E n =o(^"[/]W — fc n [f](t + h)) 2 can also be made arbitrarily small on I by taking 
|ft| sufficiently small. □ 



Lemma 4.23 ([12]). If g' 6 L 



then lim 
|ft|->o 



- 0(t + ft) 



0; thus, 



the curve g(t) = ((IC n [g}{t))) n€]N is differentiable, and (g)'(t) = ((/C"[ 5 ']0)» 



new- 



Proof. Let / be any finite interval; since g' S LjJ,, for every e > there exists 
N such that E^at+i b'l ( u ) 2 < e / 8 for a11 ^ Since functions K n [g'\ (it) 
are uniformly continuous on /, there exists a (5 > such that for all *i,*2 £ -f. if 
|*i - * 2 | < <5 then £^ =0 (/C n [g'](*i) - K. n [g'}(t 2 )) 2 < e/2. Let ft be an arbitrary 
number such that |ft| < 6; then for every t there exists a sequence of numbers 
that lie between * and t- ft, and such that {K n [g]{t)-K n [g]{t-h))/h = IC n [g'}(^ n ). 
Thus, for all * e I, 



E- 

n=0 



<?(*) - <?(* - ft) 



n 2 



</(*) 



n=0 



< 



N 

E 

n=0 



JV+1 



(^[ 5 '](^)-/C"[ 5 '](*)) 

OO 

2 £ JC n [g'](t) 2 <e/2 + 4e/8 = e 



(K n [g'}(e n )-]C n [g'}(t)]y 

OO 

2 E ^a'KO 2 



n=N+l 



□ 



Since £™[m](*) £ L 2 M for all u, if we let e k+1 (t) = (((IC k o /C n )[m](*))) neK for fc > 0, 
then by Lemma f4. 231 £?&(*) are differentiable for all fc. Since £2 is complete and ei(*) 
is continuous, e\(t) has an antiderivative H(t). Using (|26[) . we have 

ei(*) = #'(*); 

e/(*) = {{{p o £° o £")[m](f)» neM = 7 o (((/C 1 o K")[m](*))> neK 

= 7oe 2 (*); 

= -7fc- 2 (((/C fc - 2 oJC n )[m]{t))) nev + lk - 1 (((lC k o/C")[m](*))) neK 

= -7k-2 efe-i(*) + 7fc-i 4+i(*), for fc > 2. 

This means that the curve H (t) is a helix in I2 because it has constant curvatures 
Kk = 7fe— 1 for all fc > 1; the above equations are the corresponding Frenet-Serret 
formulas and efc+i(t) = (((/C fc o JC n )[m](t))) nerl for fc > form the orthonormal 
moving frame of the helix H{t). 
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5. Examples 

We now present a few examples of chromatic derivatives and chromatic expan- 
sions, associated with several classical families of orthogonal polynomials. More 
details and more examples can be found in [9]. 

5.1. Example 1: Legendre polynomials/Spherical Bessel functions. Let 

L n (oj) be the Legendre polynomials; if we set Pn(u>) = ^/2n + 1 L„(w/7r) then 

The corresponding recursion coefficients in equation (|17p are given by the for- 
mula 7„ = n(n + l)/y/A(n + l) 2 — 1; the corresponding space L 2 ,-. is L 2 [— ir,n]. 
The space LjJ, for this particular example consists of all entire functions whose 
restrictions to R belong to L 2 and which have a Fourier transform supported in 
[— 7r, 7r]. Proposition 14.191 implies that in this case our locally defined scalar prod- 
uct (f,g) M , norm \\f\\ M and convolution (/ * M g)(t) coincide with the usual scalar 
product, norm and convolution on L%. 

5.2. Example 2: Chebyshev polynomials of the first kind/Bessel func- 
tions. Let P£(to) be the family of orthonorma! polynomials obtained by normal- 
izing and rescaling the Chebyshev polynomials of the first kind, T n (ui), by setting 
P£{uj) = 1 and P£(u) = V§ T n {u/ir) for n > 0. In this case 

f F n »P^(a;) = S(n - m). 

By Proposition 14.151 the corresponding space LjJ, contains all entire functions 
fit) which have a Fourier transform f(u>) supported in [—it, tt] that also satisfies 
Vtt 2 — w 2 |/(w)| 2 du < oo. In this case the corresponding space L^, contains 
functions which do not belong to L 2 ; the corresponding function (T4T)|) is m(z) = 
Jq(ttz) and for n > 0, K. n [m](z) = (— l) n V2 J n (nz), where J„(z) is the Bessel 
function of the first kind of order n. In the recurrence relation ([26]) the coefficients 
are given by 70 = 7r/v2 and 7„ = 7r/2 for n > 0. 

The chromatic expansion of a function f(z) is the Neumann series of f{z) (see 
[23]), 

OO 

/(f) = /(u)Jo(tt(z - u)) + V2 K n [f]iu)J n (n(z - u)). 

n=l 

Thus, the chromatic expansions corresponding to various families of orthogonal 
polynomials can be seen as generalizations of the Neumann series, while the families 
of corresponding functions {A^ n [m](z)}„ e i N can be seen as generalizations and a 
uniform representation of some familiar families of special functions. 

5.3. Example 3: Hermite polynomials/Gaussian monomial functions. Let 

be the Hermite polynomials; then polynomials P"(u>) — {2 n n\)^ 1 / 2 H n {uj) 

satisfy 

/OO J 
P„»P-(o;)e- 2 ™=6in-m). 
-00 V n 
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The corresponding space LjJ, contains entire functions whose Fourier transform 
f(u>) satisfies l/( w )| 2 e ^ dw < oo. In this case the space contains non- 
bandlimited signals; the corresponding function defined by (|40|) is m{z) = e~ z / 4 
and K. n [m](z) = (— l) n (2 n n!) -1 / 2 z n e~ z / 4 . The corresponding recursion coeffi- 
cients are given by j n — W{n + 1) /2. The chromatic expansion of /(z) is just the 
Taylor expansion of f(z) e z / 4 , multiplied by e~ z / 4 . 

5.4. Example 4: Herron family. This example is a slight modification of an 
example from [9]. Let the family of orthonormal polynomials be given by the 
recursion L (uj) = 1, L 1 (u>) = ui, and L n+ i(uj) — u/(n+l)L n (u>)— n/(n+l)L n _i(w). 
Then 



1 /* /7TUJ\ 

— J L(m,oj) L(n,u>) sech^-^-J duj = S(m 



In this case m(z) — sechz and JC n [m](z) = (— l)"sechz tanh™ z. The recursion 
coefficients are given by j n = n+ 1 for all n > 0. If i?„ are the Euler numbers, then 
sechz = J2^=o z 2n /(2n)l, with the series converging only in the disc of radius 
7r/2. Thus, in this case m(z) is not an entire function. 

6. Weakly bounded moment functionals 

6.1. To study local convergence of chromatic expansions of functions which are not 
in Lj^ we found it necessary to restrict the class of chromatic moment functionals. 
The restricted class, introduced in [12], is still very broad and contains functionals 
that correspond to many classical families of orthogonal polynomials. It consists of 
functionals such that the corresponding recursion coefficients j n > appearing in 
(PS)) are such that sequences {j n }nen an d {7n+i/7ri}neiN are bounded from below 
by a positive constant, and such that the growth rate of the sequence {7n}n.ew i s 
sub-linear in n. For technical simplicity in the definition below these conditions are 
formulated using a single constant M in all of the bounds. 

Definition 6.1 ([12]). Let M be a moment functional such that for some 7„ > 
(HHJ) holds. 

(1) M. is weakly bounded if there exist some M > 1, some < p < 1 and some 
integer r > 0, such that for all n > 0, 

(59) i < 7 n <M(n + r)P, 

(60) < M 2 . 

7«+l 

(2) M is bounded if there exists some M > 1 such that for all n > 0, 

(61) Ji- ln ~ M - 

Since (|26|) is assumed to hold for some 7„ > 0, weakly bounded moment function- 
als are positive definite and symmetric. Every bounded functional M. is also weakly 
bounded with p = 0. Functionals in our Example 1 and Example 2 are bounded. 
For bounded moment functionals M. the corresponding moment distribution a(ui) 
has a finite support [2] and consequently m(t) is a band- limited signal. However, 
m(t) can be of infinite energy (i.e., not in L' 2 ) as is the case in our Example 2. 
Moment functional in Example 3 is weakly bounded but not bounded (p = 1/2); 
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the moment functional in Example 4 is not weakly bounded (p = 1). We note 
that important examples of classical orthogonal polynomials which correspond to 
weakly bounded moment functionals in fact satisfy a stronger condition from the 
following simple Lemma. 



Lemma 6.2. Let M. be such that holds for some 7„ > 0. If for some < p < 1 
the sequence j n / nP converges to a finite positive limit, then M. is weakly bounded. 

Weakly bounded moment functionals allow a useful estimation of the coefficients 
in the corresponding equations (f35|) and (f36|) relating the chromatic and the "stan- 
dard" derivatives. 

Lemma 6.3 ([12]). Assume that M. is such that for some M > l,r > and p > 
the corresponding recursion coefficients 7„ for all n satisfy inequalities (|59[) . Then 
the following inequalities hold for all k and n: 

\(K n o D k )[m](0)\ 



(62) 
(63) 





-t k ~ 






k\_ 


(0) 



< 
< 



{2M) k (k + r)\P; 
(2M) n . 



Proof. By ([30)) . it is enough to prove (|62[) for all n, k such that n < k. We proceed 
by induction on k, assuming the statement holds for all n < k. Applying ()26|1 to 
D k [m]{t) we get 



\(tC n o D k+1 )[m](t)\ < ln \(IC n+1 o D k )[m] (t)\ + ln _ 1 {(IC^ 1 o D k )[m]{t)\. 
Using the induction hypothesis and ((30)) again, we get for all n < k + 1, 

|(/C"o J D fc+1 )[m](0)| < (M(k + l + r) p + M(k + rf){2M) k {k + r)\ p 



{2M) k+ \k + l + r) 



Similarly, by (|37|) , it is enough to prove (|63|) for all k < n. This time we proceed 
by induction on n and use (|26|) . (|59|) and (|60|l to get 



/C 



< M 



n-l 



fe! 



(fc-1)!. 

By induction hypothesis and using (|37[) again, we get that for all k < n 
JC n+1 [4] (0) < M(2M) n + M 2 (2M) n - 1 < (2M) n+1 . 



1, 
□ 



Corollary 6.4 ([12]). Let M. be weakly bounded; then for every fixed n 



lim I (AC" oD k ) [m] (0)/fc! 



i/fc 



= 0, 



and t/ie convergence is uniform in n. 

Proof. Let i?(fc) = (k + r)\/k\; then i?(fc) is a polynomial of degree r, and, by (|62l . 



(64) 



(/C" o i? fc ) [m] (0) 



fc! 



i/fe 



< 



2MR{kY' k 2Me 1 ~P R(ky' k 



k\(i-p)/k 



< 



a 



Corollary 6.5 ([12]). Let m(z) correspond to a weakly bounded moment functional 
M; then 

l/k 



(65) 



lim 



(f) 



l/k 



lim 



m k {0) 
k\ 



0. 
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Thus, since (|38|) is satisfied with p = 0, every weakly bounded moment functional 
is chromatic. 



Note that this and Proposition ^. 51 imply that m(z) is an entire function. If ([59 
holds with p — 1 , then Lemma 16.31 implies only 



(66) limsup 

k — >oo 



(/C« o D k ) [m] (0) 1 



k\ 



< 2M. 



Example 4 shows that in this case the corresponding function m(z) need not be 
entire. Thus, if we are interested in chromatic expansions of entire functions, the 
upper bound in (|59[) of the definition of a weakly bounded moment functional is 
sharp. 

Lemma 16.31 and Proposition 14.71 imply the following Corollary. 

Corollary 6.6 ([12]). If M. is weakly bounded, then the corresponding family of 
polynomials {P J ^'(a;)} ne ] N is a complete system in L^^y 

Thus, we get that the Chebyshev, Legendre, Hermite and similar classical fami- 
lies of orthogonal polynomials are complete in their corresponding spaces L^f^y 

To simplify our estimates, we choose K > 1 such that for p, M and r as in 
Definition 16.11 for all k > 0, we have 



(67) {2M)Hk + rr < Rk 
K ' k\P 

The following Lemma slightly improves a result from [12]. 

Lemma 6.7. Let Ai be weakly bounded and p < 1 and K > 1 such that (|59p and 
(|67[) hold. Let also k be an integer such that k > 1/(1 — p). Then there exists a 
polynomial P(x) of degree k — 1 such that for every n and every z 6 C, 

(68) \ K n [m]{2)]< L_l_p {lz])e \K Z \ 

Proof. Using the Taylor series for K n [m]{z), (H, ((62]) and (67}, we get that for z 
such that \Kz\ > 1, 

I^HWI < E 1 lKZr n < ?E'" Zr 



< 



(n + mV. 1 ^ ~ nl 1 -? ^ to! 1 /* 
\Kz\ n ^ \ Kz \k[ m /k\+k-i 
tiP-p ^ |m/fc|! 

m=0 L ' J 



|j^|n+fc-l oo^ \Kz^_ 
n\ x -P 

If |1C2| < 1, then a similar calculation shows that for such z we have |/C™[m](z)| < 
k e \Kz\ n jri}-P. The claim now follows with P{\z\) = k^Kz^- 1 + e). □ 
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6.2. Local convergence of chromatic expansions. 

Proposition 6.8. Let M. be weakly bounded, p < 1 as in (|59[) . f(z) a function 
analytic on a domain G C C and mcG. 

(1) If the sequence |/C"(u)/n! 1_p | 1 /™ is bounded, then the chromatic expansion 
CE M [/, u](z) of f(z) converges uniformly to f{z) on a disc D C G, centered 
at u. 

(2) In particular, if \K n {u)/n\ 1 ~ p \ 1 / n converges to zero, then the chromatic ex- 
pansion CE^f/, u]{z) of f(z) converges for every z S G and the convergence 
is uniform on every finite closed disc around u, contained in G. 

Proof. Assume that R is such that limsup ?woo \K, n [f](u)/n\ 1 - p \ 1 / n < R. Then 
\rC n [f](u)\ < R n nl 1 -P for all sufficiently large n. Let K and k be such that |68|) 
holds; then Lemma l6.7l implies that for all sufficiently large n, 

\JC n [f]{u)JC n [m]{z - u)\^ n < RK(P(\z - u\)e^ K ^ u ^ k )^ n \z - u\. 

Thus, the chromatic series converges uniformly inside every disc D C G, centered 
at u, of radius less than 1/{RK). Since 

oo 

K? [CK M [f, u](z)]\ z=u = J2(-l) n K n if] («) <P ° £ n ) [m] (0) = X? [/] («) , 



CA M [/,u](z) converges to f(z) on D. 

Lemma 6.9. Let M be as in Definition \6.1{ Then 



□ 



lim sup 



X n [f](t 



l/n 



< 2M lim sup 



Proof. Let /3 > be any number such that limsup,^^ \f^(u)/n\ l p \ x / n < /3; 
then there exists Bp > 1 such that \f {k) {u)\ < Bp k\ l - p (3 k for all k. Using 
and (|6"5|) we get 



\fC n [f](u)\ 



< 



k=Q 



K, n 



(0) 



|/ (fe) NI < (2Af)"S /3 y" fc! 1 "^' 



< (2MrB f3 j2 - 



/3 



fe=0 



Summation of the last series shows that \K n [f](u)\ < 2B p (2M(3) n n\ 1 -"P for suffi- 
ciently large n. □ 



Corollary 6.10. Let M be weakly bounded, p < 1 as in (|59[) . f(z) a function 
analytic on a domain G C C and u 6 G. 

fij sequence |/^ n - ) (u)/n! 1 ~ p | 1 / Tl is bounded, then the chromatic expansion 

CET[/,«](z) of f{z) converges uniformly to f(z) on a disc D C G centered 

at u . 

(2) In particular, if |/(")(u)/n! 1_p | 1 /™ converges to zero, then the chromatic 
expansion CE^/, u](z) of f(z) converges for all z G G and the convergence 
is uniform on every closed disc around u, contained in G. 

Corollary 6.11 ([12]). If M. is bounded, then for every entire function f and all 
u, z S C. the chromatic expansion CE[f, u](z) converges to f(z) for all z, and the 
convergence is uniform on every disc around u of finite radius. 
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Proof. If f(z) is entire, then for every u, limn-™ = 0. The Corollary 

now follows from Corollary 16. 101 with p = 0. □ 

Proposition 6.12. Assume M. is weakly bounded and let < p < 1 be such that 
(|59p holds, and k such that k > 1/(1 — jj). TTien there exists C, L > smc/i i/iai 
!/(*)! <C||/|Le^l* /or flll/(z)£L J 



Proof. Since f(z)E LjJ,, the chromatic expansion of /(z) and (|68|) yield 

/ oo oo \ 

< ^|/C"[/](0)| 2 ^|/C"[m](z)| 2 



< ll/IL^(M) e 



A'.; 




which, using the method from the proof of Lemma [6.71 can easily be shown to imply 
our claim. □ 

Note that for bounded moment functionals, such as those corresponding to the 
Legendre or the Chebyshev polynomials, we have p = 0; thus, Proposition 16.1 21 im- 
plies that functions which are in L^, are of exponential type. For M. corresponding 
to the Hermite polynomials p=l/2 (see Example 3); thus, we get that there exists 
C.L > such that \f(z)\ < C \\f\\ M e L \ z \ 2 for all / e \? M . It would be interesting 
to establish when the reverse implication is true and thus obtain a generalization 
of the Paley- Wiener Theorem for functions satisfying \f(z)\ < Ce L l z l for k > 1. 

6.3. Generalizations of some classical equalities for the Bessel functions. 

Corollaries 16 . 1 1 1 and 16.101 generalize the classic result that every entire function can 
be expressed as a Neumann series of Bessel functions [23], by replacing the Neu- 
mann series with a chromatic expansion that corresponds to any (weakly) bounded 
moment functional. Thus, many classical results on Bessel functions from [23] 
immediately follow from Corollary I6.11[ and, using Corollary \6.10\ generalize to 
functions JC n [m](z) corresponding to any weakly bounded moment functional M. 
Below we give a few illustrative examples. 

Corollary 6.13. Let P^{lj) be the orthonormal polynomials associated with a 
weakly bounded moment functional M; then for every z£C, 



(69) e'-^^HrHiz). 

Proof. If p < 1 then 



n=0 



liin — — - — = lim — LJ— = 0, 

n p n^oo n p 



and the claim follows from Proposition 16.101 and (f27|) . □ 

Corollarv l6 . 131 generalizes the well known equality for the Chebyshev polynomials 
T„(u) and the Bessel functions Jn(z), i.e., 

oo 

c^ z = J (z) + 2 i n T n (w)J n {z). 
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In Example 3. 16.131 becomes the equality for the Hermite polynomials H n (u>): 

H n (uj) ( iz s 



3 iwz __ ' 



n=0 



1)1 



! V 2 



By applying Corollary 16.101 to the constant function f(t) = 1, we get that its 
chromatic expansion yields that for all z 6 C 

oo / n \ 

with 7 n the recursion coefficients from ()17p . This equality generalizes the equality 

OO 

J (z) + 2^ J 2n (z) = 1. 

Using Proposition 14.161 to expand m(z + u) £ LjJ, into chromatic series around 
z = 0, we get that for all z, it 6 C 

oo 

m(z + w) = ^(-l) ,l /C n [m]( U )/C"[m](z), 

n=0 

which generalizes the equality 

oo 

J (z + «) = Jo(w)Jo(^) + 2 ^(-l) n J„(u)J„(z). 



7. Some Non-Separable Spaces 

7.1. Let .M be weakly bounded; then periodic functions do not belong to be- 
cause X)^o ^ n if] (^) 2 diverges. We now introduce some nonseparable inner product 
spaces in which pure harmonic oscillations have finite norm and are pairwise or- 
thogonal. 

Note: In the remainder of this paper we consider only weakly bounded moment 
junctionals and real functions which are restrictions of entire functions. 



Definition 7.1. Let let < p < 1 be as in (|59|) . We denote by C M the vector space 
of functions such that the sequence 

1 n 

(70) 4(t) = 7 — nl - P E^[/]W 2 

1 ' k=0 

converges uniformly on every finite interval ICR. 
Proposition 7.2. Let f,geC M and 

n 

(71) = fw+1 u- P £^ fc [/](*)^b](*); 

then the sequence {c^ 9 (i)}ne]N converges to a constant function. In particular, 
{ v L{ty\n£m a ^ so converges to a constant function. 
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Proof. Since v^it) and v^it) given by (|70|) converge uniformly on every finite in- 
terval, the same holds for the sequence <r^ 9 (i). Consequently, it is enough to show 
that for all t, the derivative ^(^(t) 1 of ^(^(t) satisfies linv^oo cr^ 9 (ty = 0. Let 

S k (t) = K, k [f](t) 2 + JC k+1 [f}(t) 2 + K k [g]{tf +K k+1 [g]{tf- 

then, since f,g£ C M , the sequence l/(n + J2k=o &k(t) converges everywhere 
to some a(t). We now show that if t is such that a(t) > 0, then there are infinitely 
many k such that Sk{t) < 2a(t)k~ p . Assume opposite, and let K be such that 
S k (t) > 2a(t)k~P for all k > K. Then, since 

A r +1 __, (n + iy-p-x 1 -? 

k p > / x p dx = — 



k=K 



K 



1-p 



we would have that for all n > K , 



(n+iy-P ~ (n + iy-P (n + -p) 

However, since < p < 1, this would imply Y^k=o ^k(t)/( n + > for 

all sufficiently large n, which contradicts the definition of a(t). Consequently, for 
infinitely many n all four summands in S n (t) must be smaller than 2a(t) n~ p . For 
those values of n we have 

\JC n+1 [f]{t)K n [g]{t)\ + \IC n [f](t)JC n+1 [g}(t)\ < 4a(t)n^. 

Since M. is weakly bounded, (|3ip and (|59p imply that for some M > 1 and an 
integer r, 

H 9 (t)'\ < ^^ (\)C n+1 [f}(t)}C n [9}(t)\ + |/C"[/]ft)/C" +1 [g](t)|). 

Thus, for infinitely many n we have 

, f , 4M(n + r)P n'* a(t) 
' ' (n + l) L p 

Consequently, liminfn^oo |cr,;{ 9 (t)'| = and since lim^^oo o^ g (t)' exists, it must be 
equal to zero. □ 

Corollary 7.3. Let C^ 1 be the vector space consisting of functions f(t) such that 
limn—xx, (t) — 0; then in the quotient space = C M /C^ we can introduce a 
scalar product by the following formula whose right hand side is independent oft: 

1 - 

(72) (f,g)» = lim - V/C fc [/](t)/C fc [ g ]ft). 

n— >oo n + 1 1 P z — ' 

V ' k=0 

The corresponding norm on C^ 4 is denoted by || ■ || . Clearly, all real valued 
functions from belong to Cq 4 . 

Proposition 7.4. If f € C^ 4 , then the chromatic expansion of f{t) converges to 
f(t) for every t and the convergence is uniform on every finite interval. 
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Proof. Since l/(n + l) 1 ^ Efc=o ^ fc [/](*) 2 converges to < (||/||") a < oo, for all 
sufficiently large n, 

/ n \ V(2n) 

l/c-i/iwi 1 /" < £> fe [/](*) 2 

\fe=o / 
< (2||/|r) 1/n (n + f) (1 ~ p)/(2 ™ ) . 

Thus, |/C n [/](i)/n! 1 - p | 1/ ™ -> 0, and the claim follows from Proposition EH □ 

M 

Since K. n [m](t) e L^,, we have E^oC-^^'I/K ) ^ J 'H(*) = for all n. 
Thus, the chromatic expansion of / € C^ 1 does not converge to /(i) in C2 4 . More- 
over, there can be no such series representation of functions / G , converging in 
Cjj 4 , because the space C^ 1 is in general nonseparable, as the remaining part of this 
paper shows. 

7.2. Space C^ 1 associated with the Chebyshev polynomials (Example 2). 

For this case the corresponding space C^ 1 will be denoted by CJ , and in (|59|) we 
have p = 0. Thus, the scalar product on CJ is defined by 



n 

(/,«,)*■= hm ^/C fe [/](t)/C fe [ 5 ](t). 

n— >oo n + 1 * — ' 



fe=0 



Proposition 7.5. Functions f u {t) = \[2 sin and <? w (t) = \J~2 cos for < lu < 

7r /orm an orthonormal system of vectors in . 

Proof. From (|2T|) we get 

Efc=0 P 2fe M P 2fc <» sin ^ sin CT * 



2n + l 



^ + Um ELp P 2 T fc+l Mjafc+l (gj COS Uft COS trt 

n— >oo 2n + 1 

Since -P^(w) < \/2 on (0, 7r), fTS]) implies that for w, cr g (0, 7r) and uj 7^ cr, 

lim ELg^MW = lim = 

n— >oo n + 1 n— >oo 2(n + f ) 

Since PJ n (uj) are even functions and -P^+i^) odd, this also implies that 

lim ELo^H^W = Um EL"c^ + i(^ + iW _ 

n— >oo 2n + f n— >oo 2n + 1 

Thus, by J73]), (f UJ ,f a } M = 0. Using (J2D]), one can verify that for < w < tt 

I V^pT/ \a 1 + 2n sin((2n + 1) arccosw) _^ 
n+ f ^ k ^> ~ 2n + 2 + (2n + 2)VT^Z? " 

Thus, ([19]) implies that for < w < tt 

n n 

lim ^— tE^h 2= lim ^— tE^+iH 2 = o 

7woo 2n + 1 ' — ' n— >oo 2n + 1 * — ' 2 

fc=0 fe=0 

I I / — 1 1 M 

Consequently, v2sinu;£ =1. □ 
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7.3. Space Cf 4 associated with the Hermite polynomials (Example 3). The 

corresponding space in this case is denoted by Cf , and in |59|) we have p = 1/2. 
Thus, the scalar product in Cf is defined by 

n 

Proposition 7.6. For all lu > functions fu(t) — sinujt and <? w (t) = cosut form 
an orthogonal system m Cf , and \\fu\\ M = \\9u\\ M — c " ' 2 /v2tt- 

Proof. For all w and for n — * oo, 

„, , r(n + l) 5 ^ 2 / , rwr\ 

Rf (w) „ v , ; - e— cos ( V2n + 1 w ) -»• 0; 

" V 7 2fr(f + 1) V v 2 7' 

see, for example, 8.22.8 in [18]. Using the Stirling formula we get 

i 

- e— cos (^v^n + 1 w - —j -> 0. 

This fact and lfl8|) are easily seen to imply that (fu,fcr) M — for all distinct 
ui,a > 0, while $73J) and ((T5J) imply that 

2 P a r Cl) )2\ 



2 fe=0 ^2fc ( w ) 2 S fc=0 ^2fc+ 1 ( 



(75) lim ^ fc = u 2fcv y - ^"=" " K+1 V y - = 0. 

1 V2n + 1 V2ri + 1 / 



Since ff^(w) = 2 ri H n _i(w), we have P£(u)' = V^n P%_ x (u). Using this, (J7U) 
and (|2"0|) one can verify that 



which, together with (75]), implies that ||/ ta ,||' M = / 2 / \/2tt. □ 

Note that in this case, unlike the case of the family associated with the Chebyshev 
polynomials, the norm of a pure harmonic oscillation of unit amplitude depends on 
its frequency. 

One can verify that propositions similar to Proposition 17.51 and Proposition 17.61 
hold for other classical families of orthogonal polynomials, such as the Legendrc 
polynomials. Our numerical tests indicate that the following conjecture is trueQ 

Conjecture 7.7. Assume that for some p < 1 the recursion coefficients j n in (fT7|) 
are such that "f n /n p converges to a finite positive limit. Then, for the corresponding 
family of orthogonal polynomials we have 

1 " 

(76) < lim - P " < 00 

n— >oo (n + 1 1 p L — ' 

for all to in the support sp(a) of the corresponding moment distribution function 
a(u>). Thus, in the corresponding space Cf 4 all pure harmonic oscillations with 
positive frequencies u> belonging to the support of the moment distribution a(u>) 
have finite positive norm and are mutually orthogonal. 



'We have tested this Conjecture numerically, by setting 7„ = n p for several values of p < 1, and 
in all cases a finite limit appeared to exist. Paul Nevai has informed us that the special case of 
this Conjecture for p = is known as Nevai-Totik Conjecture, and is still open. 
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Note that (f20|) implies that (Tf6|) is equivalent to 

< hm -, . , „ < oo. 

n^oo (n + l) 1 - 2 ^ 

8. Remarks 

The special case of the chromatic derivatives presented in Example 2 were first 
introduced in [11]; the corresponding chromatic expansions were subsequently in- 
troduced in [14]. These concepts emerged in the course of the author's design of a 
pulse width modulation power amplifier. The research team of the author's startup, 
Kromos Technology Inc., extended these notions to various systems corresponding 
to several classical families of orthogonal polynomials [5,9]. We also designed and 
implemented a channel equalizer [8] and a digital transceiver (unpublished), based 
on chromatic expansions. A novel image compression method motivated by chro- 
matic expansions was developed in [3,4]. In [6] chromatic expansions were related 
to the work of Papoulis [17] and Vaidyanathan [19]. In [15] and [20] the theory was 
cast in the framework commonly used in signal processing. Chromatic expansions 
were also studied in [5], [1] and and [22]. Local convergence of chromatic expan- 
sions was studied in [12]; local approximations based on trigonometric functions 
were introduced in [13]. A generalization of chromatic derivatives, with the prolate 
spheroidal wave functions replacing orthogonal polynomials, was introduced in [21]; 
the theory was also extended to the prolate spheroidal wavelet series that combine 
chromatic series with sampling series. 



Note: Some Kromos technical reports and some manuscripts can be found at the 
author's web site http://www.cse.unsw.edu.au/~ignjat/diff/. 
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